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We report on progress in our study of high temperature QCD with three flavors of improved staggered quarks. 
Simulations are being carried out with three degenerate quarks with masses less than or equal to the strange 
quark mass, m s , and with degenerate up and down quarks with masses in the range 0.1m, < m u ,d < 0.6m s , 
and the strange quark mass fixed near its physical value. For the quark masses studied to date we find rapid 
crossovers, which sharpen as the quark mass is reduced, rather than bona fide phase transitions. 



The MILC Collaboration is studying high tem- 
perature QCD with three flavors of improved 
staggered quarks [T] using a one-loop Symanzik 
improved gauge action 2 and the Asqtad quark 
action 3 . This action is particularly well suited 
for the study of high temperature QCD because 
it has excellent scaling properties in the lattice 
spacing, significantly better dispersion relations 
for quarks and hadrons than standard actions, 
and substantially smaller taste symmetry viola- 
tions than the Kogut-Susskind action. 

We are considering two cases: 1) all three 
quarks have the same mass, m q ; and 2) the two 
lightest quarks have equal mass, m u ^, and the 
mass of the third quark is fixed at that of the 
strange quark, m s . We refer to these cases as 
Nf — 3 and Nf = 2 + 1, respectively. This year 
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we have reduced the quark mass in the Nf = 3 
study to m q — 0.2 m s , and that in the Nf = 2 + 1 
study to m u ,d = 0.1 m s . For the masses we have 
studied to date, we find rapid crossovers, which 
sharpen as the quark mass is reduced, rather than 
bona fide phase transitions. 

We would like to keep the physical quark 
masses fixed as we vary the lattice spacing, and 
therefore the temperature. To this end, for the 
Nf = 3 case we work along curves of con- 
stant pseudoscalar to vector mass ratio, m^/m^. 
These curves are determined from interpolations 
between the results of spectrum calculations at 
lattices spacings 0.12 and 0.20 fm. For Nf = 2 + 1 
we make use of spectrum calculations at lattice 
spacings 0.09, 0.12 and 0.20 fm to determine 
curves of fixed m^/mp and m nss /m^. Here rj ss 
and 4> are the pseudoscalar and vector mesons 
made up solely of strange quarks. The lattice 
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spacing, and hence the temperature, are deter- 
mined from the heavy quark potential for values 
of the gauge coupling and quark masses for which 
we have performed spectrum calculations, and by 
interpolation of the lattice spacing for other val- 
ues of the gauge coupling and quark mass. 

For three equal mass quarks, Nf = 3, we have 
carried out thermodynamics studies on lattices 
with four, six and eight times slices, and aspect 
ratio N s /N t = 2. Here N s and N t are the spatial 
and temporal dimensions of the lattice in units 
of the lattice spacing. To date we have studied 
quark masses in the range 0.2m s < m q < m s . 
In Fig. ^ w e plot the chiral order parameter, 
(ipip) as a function of temperature on 16 3 x 8 
lattices. The bursts are linear extrapolations of 
langleiptp) from the runs with the two lightest 
quark masses to m q = for fixed temperature. 
This figure suggests that for m q = there is un- 
likely to be a phase transition for temperatures 
above 175 MeV, but one could occur at or below 
that value. For the quark masses we have stud- 
ied, the data indicate a crossover, which sharpens 
as the quark mass decreases. The absence of a 
phase transition at these quark masses is consis- 
tent with the findings of the Bielefeld group 0]. 
The sharpening of the crossover with decreasing 
quark mass is seen more clearly in Fig. where 
we plot the total iptp susceptibility as a function 
of temperature on 12 3 x 6 lattices. 

The Nf = 2 + 1 thermodynamics studies were 
carried out primarily on 12 3 x 6 and 16 3 x 8 lat- 
tices. In addition, a number of runs were made 
on 18 3 x 6 lattices to check for finite size effects, 
which turned out to be negligible. In this phase 
of the work we performed simulations with two 
degenerate light quarks in the range 0.1 m s < 
m,u,d < 0.6 m s . In Fig. [3] we show the chiral or- 
der parameter as a function of temperature on 
16 3 x 8 lattices. Here the bursts are a linear ex- 
trapolation from the runs with the two smallest 
light quark masses to m Ui d = for fixed tempera- 
ture and heavy quark mass. This figure indicates 
that there is unlikely to be a phase transition for 
mu,d = and the heavy quark mass equal to that 
of the strange quark for T > 175 MeV. The sharp- 
ening of the crossover as m Ut d decreases is clear. 

Quark number susceptibilities |5I6| are of par- 
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Figure 1. The chiral order parameter, {ipij)), on 
16 3 x 8 lattices for Nf — 3. The bursts are lin- 
ear extrapolations of (tptp) to m q — at fixed 
temperature. 
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Figure 2. The t/jtf> susceptibility as a function of 
temperature on 12 3 x 6 lattices for Nf = 3. 
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Figure 3. The chiral order parameter, {ipip), on 
16 3 x 8 lattices for Nf = 2 + 1. The bursts are 
linear extrapolations of (ipip) to m q = at fixed 
temperature. 
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Figure 4. The triplet and strange quark number 
susceptibilities as a function of temperature for 
two light quarks with mass 0.2 m s and one heavy 
quark with mass m s . Results are shown for 12 3 x6 
and 16 3 x 8 lattices. The lines on the right of the 
figure indicate the values for free quarks in the 
continuum, and on the finite lattices on which 
the simulations were carried out. 



ticular interest because they are related to event 
by event fluctuations in heavy ion collisions, and 
can therefore be directly measured |Zj . They also 
provide excellent signals for the crossover. In Fig- 
ure 01 we show results for the triplet and strange 
quark number susceptibilities on both 12 3 x 6 and 
16 3 x 8 lattices. The close agreement between the 
Nt = 6 and 8 results illustrates the excellent scal- 
ing properties of the action in the lattice spacing. 

Of particular importance for the coming year 
will be the determination of the end point of the 
line of first order phase transitions expected for 
the Nf — 3 case, and the completion of the Nf = 
2 + 1 study at m u ,d/m s = 0.1. 
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